Interplane magnetic coupling effects in the multilattice compound Y 2 Ba 4 Cu 7 Oi5 
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We investigate the interplane magnetic coupling of the multilattice compound Y2Ba4Cu70i5by 
means of a bilayer Hubbard model with inequivalent planes. We evaluate the spin response, effective 
interaction and the intra- and interplane spin-spin relaxation times within the fluctuation exchange 
approximation. We show that strong in-plane antiferromagnetic fluctuations are responsible for a 
magnetic coupling between the planes, which in turns leads to a tendency of the fluctuation in the 
two planes to equalize. This equalization effect grows whit increasing in-plane antiferromagnetic 
fluctuations, i. e., with decreasing temperature and decreasing doping, while it is completely absent 
when the in-layer correlation length becomes of the order of one lattice spacing. Our results provide 
a good qualitative description of NMR and NQR experiments in Y2Ba4Cu70is. 
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Although many models for high-Tc cuprates are re- 
stricted to a single layer, it has become clear that both 
superconducting and magnetic properties of these mate- 
rials are affected by the coupling between two or more 
layers. A rather strong coupling between the layers has 
been observed principally by inelastic neutxcoi scatter- 
ings (INS) and nuclear magnetic resonanceoO (NMR). 
Furthermore, the observation of a qualitatively different 
behavior of the odd and even channel in INScI and of a 
bilayer splitting of the Fermi surface found in- angular 
resolved photoemission experiments (ARPES)BeI demon- 
strate that low energy excitations of cuprates are affected 
by the presence of more than one layer per unit cell. An 
exciting perspective on the nature of the coupling be- 
tween CuOvlayers was offered by NMR experiments by 
Stern et al. on Y2Ba4Cu70is (247). This material has a 
variety of structural similarities to the extensively studied 
YBa 2 Cu 3 7 (123) and YBa 2 Cu 4 8 (124) systems. The 
compound 247 can be considered as a natural multilat- 
tice, whose bilayers are build up of one Cu02 layer which 
belongs to the 123 block and one layer to the 124 block. 
Based on the analysis of the NQR spectra it turned out 
that the charge carrier content in these nonequivalent 
adjacent layers is very close to that of the related parent 
compounds of the two blocks, 123 and 247. Interestingly, 
the highest transition temperature (T c = 95 K) occurs in 
the 247 compound, in comparison with the 92 K of 123 
and 82 K of the 124 system. 

In this paper, we want to provide a theoretical under- 
standing in terms of a microscopic model of some strik- 
ing experimental observations of Refs. |^,^, namely: (i) 
the spin-spin relaxation rates T^ of the two layers in 
Y2Ba4Cu70i5, measured in a spin-echo double resonance 
experiment, behave very similarly as a function of tem- 
perature, despite the different doping of the layers; (ii) 
the spin-spin relaxation rate in the 124 (247) layer of 
Y2Ba4Cu70i5is reduced (enhanced) with respect to one 
of the constituent compound at low temperatures; (iii) 
the interplane transverse relaxation rate, increases for 
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decreasing temperature faster than the intraplane one. 
The overall features that can be inferred from these ex- 
periments are that, for high temperatures the magnetic 
fluctuations in the two layers of Y 2 Ba4Cu 7 0i5are dis- 
connected and each layer behaves similarly to the corre- 
sponding parent compounds, whereas for decreasing tem- 
peratures, the increasing interlayer magnetic coupling en- 
forces even the slightly overdoped plane to behave like an 
underdoped system and vice-versa. 

To describe the strong electronic correlations in the 
system Y2Ba4Cu70i5, consisting of two layers with dif- 
ferent charge carrier concentration, we extend the stan- 
dard single-band Hubbard model to include two layers 
coupled by an interplane hopping matrix element t±. 
Furthermore, to produce a different hole concentration in 
the planes, we introduce an on-site energy 8 in the second 
plane. After Fourier transformation of the kinetic part 
the Hamiltonian reads: 
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where c\. z a creates a particle with spin a and mo- 
mentum k in layer l\, and i runs over the sites of the 
two layers. In contrast to the standard notation in 
the monolayer case, ek is replaced by the 2x2 ma- 
trix. H Q (k). whose components are [_ff (k)]n = — fj,, 
[H (k)} 22 =e k + 5-fi and [H (k)} 12 = [H (k)} 21 = t±. 
The bare dispersion £k includes second (t') and third (t") 
nearest-neighbor hopping processes to tetter model the 
Fermi surface of the cuprates (see, e.g.ErEj). The prop- 
erties of the interacting system are deduced from the 
Green's function and form the dynamic two-particle sus- 
ceptibilities in the framework of the fluctuation exchange 
approximation (FLEX)Ej. In this approximation, bub- 
ble and ladder diagrams are summed up in infinite or- 
der and the resulting coupled set of equations is solved 
self-consistently. Although this approximation is suit- 
able to describe magnetic fluctuations and band shadow 
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featured it fails to reproduce the pseudogap behavior in 
the spectral function and in the magnetic excitations at 
low temperatures and low dopingpil, probably due to the 
lack of conservation at the two-particle level and to the 
neglection of interference effects between particle-particle 
and particle-hole channel. Moreover, on the quantita- 
tive level, antiferromagnetic fluctuations appear to be 
underestimatecO with respect to exact results. 

In order to understand the systematics and the pa- 
rameter dependence of the intcrplanc magnetic coupling, 
we first focus our attention on two systems with differ- 
ent parameter sets t'/t, t" jt corresponding to strong and 
weak antiferromagnetic fluctuations. Specifically, we in- 
troduce a first system, for simplicity labeled by "A" , with 
t' = -0.38t,t" = -0.06i, , and a second one ("B") 
with t' = -0.20<, t' = 0.15t. System A shows much 
stronger antiferromagnetic fluctuations due to the un- 
derlying Fermi surface in comparison to B. The Hubbard 
interaction takes an intermediate value U — 4i, as ap- 
projaria.te-.for a perturbative calculation, and t± = OAt 
(cf .Eall3o) . Within the self-consistency cycle, we fix the 
on-site energy 6 and the particle number n\ = 1 — x\ 
of the first layer, while the chemical potential [i and the 
particle number of the second plane = 1 — X2 are 
adjusted at each step. Keeping this way the doping of 
the first plane fixed and changing only the doping of the 
second one, we investigate a possible magnetic coupling 
between the layers. Our question is whether the mag- 
netic fluctuations of the first plane are influenced by the 
doping of the second plane and vice-versa. 

With this in mind, we consider in Fig. [l] the static 
spin-spin function \ii (q ; w = 0) along the standard path 
(0, 0) -> (tt, 0) -> (tt, tt) -> (0, 0) in the Brillouin zone for 
the two layers (I = 1 and I — 2) and both systems A and 
B. 
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FIG. 1. Static spin susceptibility XiiiQi w 



(Jl,7t) 



(0,0) 



0) along the 

standard path in the Brillouin zone for the two layers and both 
systems A and B (see text). (The temperature is T = 0.02t, 
and xi = 0.08) 

As one can clearly see from Fig. [j](a) and (b) for the 
case of system A, the spin response in the two planes 
is strongly peaked at the antiferromagnetic wave vector 



Q = (7r, tt) indicating pronounced antiferromagnetic fluc- 
tuations in the Hubbard planes. The striking result is 
that data for a plane with a given doping (li = 1) also 
depend on whether this plane is coupled with an equiv- 
alent one or with a more or less doped one, depending 
of the value of <5. This behavior clearly suggests a strong 
magnetic connection between the planes. The decrease 
of the spin susceptibility with increasing 5 [see Fig. |l](a) 
and (b)] is due to the increased hole concentration of 
the total system moving it further away from half filling 
where antiferromagnetism is most pronounced. The sit- 
uation is rather different for system B, whose data are 
shown in Fig. |l](c) and (d). Here, a variation of S influ- 
ences the susceptibility of the second plane only, whereas 
the first plane, with constant charge carrier concentra- 
tion, is almost not affected at all. Thus, the magnetic 
fluctuations in the first plane are disconnected from the 
fluctuations in the second one. Only for a system with 
sufficiently strong antiferromagnetic correlations and a 
correlation length larger than a few lattice constants, the 
two inequivalent layers turn out to be strongly coupled, 
as in case A. 
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FIG. 2. Effective interaction 14fr(Q,w = 0) vs. T for a 
system with inequivalent layers {x\ — 0.08 and X2 = 0.11) 
in comparison with the corresponding bilayer systems with 
x\ = xi = 0.08 and xi = Xi = 0.11. The upper panel shows 
the results for system A and the lower for B. 

We now ask to what extent each layer of a system 
with inequivalent layers compares to the corresponding 
system with equivalent layers (i.e., 6 = 0). Therefore, 
we show in Fig. || the effective interaction V e g (Q, ui = 0) 
at the antiferromagnetic wave vector Q = (tt, tt) as a 
function of temperature for systems with inequivalent 
layers (xi = 0.08, x-i = 0.11,(5 = 0.2) in comparison 
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with the two corresponding systems with identical lay- 
ers (xi = X2 = 0.08) and (xi = X2 = 0.11). In the upper 
panel of Fig. || we show the results for system A. Here, we 
clearly observe considerable differences between each of 
the layers in the inequivalent-layer system with respect 
to its counterparts in the equivalent-layer system at low 
temperatures. Specifically, the effective interactions in 
the two layers of the system with inequivalent layers are 
"enclosed" between the ones in the corresponding layers 
of the system with identical layers. Thus, the interlayer 
magnetic correlations lead to a tendency for the mag- 
netic fluctuations of the two planes to equalize although 
their doping is different. The results for system B are 
strikingly different as can be seen in the lower panel of 
Fig. ||. Here, the two planes with different dopings in the 
inequivalent-layer system behave essentially like the cor- 
respondingly doped layers in the system with equivalent 
layers. For this less magnetic system the fluctuations in 
the two layers are effectively decoupled and do not in- 
fluence each other. The results displayed in Fig. || thus 
again show that the antiferromagnetic correlations within 
the layers lead to a partial equalization of the magnetic 
fluctuations in the two planes. 

A similar magnetic equalization effect has been ob- 
served by Scalettar et al.tZl in Quantum-Monte-Carlo 
simulations of an half-filled layer coupled to a doped one. 
According to these authors, this equalization effect can 
be explained by perturbation in t± as a competition be- 
tween virtual hopping processes between the two layers 
with energy scale <xt±/5 and exchange energy J oc t 2 /U. 
We believe that our use of the FLEX approximation is 
essential in order to additionally explain the crossover 
from the "equalized" regime (A) to the disconnected one 
(B), which accounts for the qualitative behavior of the 
experiments on Y2Ba4Cu70i5U as we will show below. 
This is due to the feedback effect of the Green's func- 
tions renormalization in the FLEX approximation. 

The equalization effect observed in system A qual- 
itatively describes the experimental situation for the 
inequivalent-layer system YJj^CuyOis. Indeed, NMR 
experiments by Stern et al.fj measuring I/T2G for the 
two layers of Y2Ba4Cri70i5 and for the two correspond- 
ing systems with equivalent layers, YBa 2 Cu307 and 
YBa2Cu40g, show a similar equalization tendency as in 
Fig. @. To further demonstrate the similarities of our cal- 
culations with the experimental results, we also evaluate 
the nuclear spin-spin relaxation time T2G for spins within 
the two planes. The generalized relaxation rate T^ G ~ 1 
measuring the interaction between a spin in plane I and 
one in V is related to a weighted sum over the static spin 
susceptibility X11 (Qi w ~ PiHfr 3 ^ re * ne momenta q = Q 
have the strongest weight BEjEJ. In our calculations, we 
assume that the hyperfine coupling constants entering 
the form factor are the same for both-olanes and use the 
values given by Barzykin and PinesEil. Hence possible 
difference in the relaxation times of the two planes can 
solely arise from different spin responses in the layers. 



In a previous work on Y2Ba4CuyOi5, Millis and Monien, 
who determined the value of the interlayer exchange in- 
teraction from the experimental data of Reffl, have con- 
sidered an alternative point of view and neglected the dif- 
ference between the spin susceptibilities of the two layeis 
by considering different hyperfine coupling constants^. 
This approach is justified by our results which show for 
systems with strong antiferromagnetic correlations a pro- 
nounced coupling of the magnetic response of the two 
layers and thus a tendency to equalization of the suscep- 
tibilities. This primarily causes a similar temperature 
dependence of the two in plane susceptibilities. Never- 
theless, quantitatively, the magnitudes of the magnetic 
response of the two layers stays different. 
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FIG. 3. Temperature dependence of the Gaussian com- 
ponent of the spin-spin relaxation rate T2G of a bilayer sys- 
tem consisting of two different layers with doping x\ — 0.16 
and X2 = 0.24, respectively (S = 0.4). For comparison, we 
show also the result for two bilayer systems with equivalent 
layers with doping x\ — X2 — 0.16 and xi = X2 = 0.24, 
respectively. The other parameters for all curves are: 
t' = -0.38t,t" = -0.06*. The inset shows the ratio T$ a /T% G 
for the case of inequivalent layers. 

In Fig. we show the in-plane I/T2G as a function 
of temperature for both planes of the system with in- 
equivalent layers again in comparison with the data for 
two corresponding bilayer systems with equivalent lay- 
ers. In this figure, the filled (open) symbols represent 
the system with inequivalent (equivalent) layers. Fur- 
thermore, the squares are related to the data for planes 
with a doping x = 0.16 and the bullets planes with dop- 
ing x = 0.24. Since the main contribution to I/J2G stems 
from x zz (Q,0), a system with rather strong magnetism 
exhibits large relaxation rates I/T^g- This explains why 
the layers with x — 0.16 show a smaller relaxation rate 
compared to x — 0.24. As expected form the results of 
Fig. ||, and as observed in NMR experiments at lower 
temperaturesB, the heavily doped plane (x% — 0.24) of 
the system with inequivalent layers shows stronger mag- 
netic fluctuations than the plane in the corresponding 
system with equivalent layers (xi = X2 = 0.24). Simi- 
larly, the magnetism of the lower doped plane (x\ = 0.16) 
is reduced with respect to the corresponding equivalent- 
layer system [x\ — X2 = 0.16) due to the coupling to 
a more doped plane. As already observed in the results 



3 



for the effective interaction shown above, the magnetic 
fluctuations of the two inequivalent planes with differ- 
ent carrier concentration tend to be partially equalized 
by interplane coupling effects. The theoretical results 
in Fig. ^, as well as the experimental results, show a 
strong increase of I/T2G with decreasing temperature. 
However, the experimentally observed decrease of I/T2G 
below T sg « 1 00K due to the opening of a pseudogap in 
the spin excitation spectrum, is not reproduced by the 
FLEX approximation, as discussed above. 

As NMR experiments show, the spin-lattice relax- 
ation rate I/T2G has the same temperature dependence 
in the two planes of Y2Ba4Cu7Oi50, the ratio R — 
(l/T 2 1 g 4 )/(f/T 2 1 ( ?, 3 ) being temperature independent and 
approximately R w 1.4—1.5. This ratio corresponds to 
R = T^q/T^q in our work, since the 124 layer in the 
coupled layer structure of Y 2 Ba4Cu70is is the one with 
lower doping (here labeled by "1") and the CuCVlayer 
from the 123 block corresponds to the second plane in our 
theoretical study. The theoretical values for R are shown 
in the inset of Fig. | as a function of temperature. For 
the parameter set chosen (corresponding to system A), 
R is approximately 1.2 — 1.4 and almost temperature- 
independent, in good agreement with the experimental 
finding. 
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FIG. 4. Intra- and inter-plane spin-spin relaxation 
rates I/T2G an d V^2Gj respectively, vs. T. The in- 
set shows the ratio between the two relaxation rates 
(asi = 0.16, x 2 = 0.24,(5 = 0.4). 

In order to have a measure of the interplane spin 
coupling in Y 2 Ba4Cu70i5, Stern et alu have carried 
out NQR spin-echo double-resonance (SEDOR) measure- 
ments, as suggested by Monien and RiceO. These experi- 
ments allow the determination of the inter-layer spin-spin 
relaxation time T\ G . The apparent feature in their results 
is that the inter-plane relaxation rate, although smaller 
than the in-plane one, increases faster for decreasing tem- 
perature, as seen from the temperature dependence of the 
ratio i?sEDOR.(T) = T^q/T^q. For the sake of compar- 
ison, we present our theoretical results for the in-plane 
relaxation rate I/T^q together with the inter-plane one 
1/To^ in Fig. |], the parameters being the same as in 
Fig. g. The fact that 1/T$ grows faster with decreasing 
temperature is even more apparent in the inset. Here, we 



show the temperature dependence of the ratio between 
the two relaxation rates. Our theoretical results thus 
reproduce the qualitative behavior observed experimen- 
tally. However, for T smaller than the spin gap the exper- 
iments show a saturation effect, which obviously cannot 
be reproduced within our approximation. 

In summary, we have studied the magnetic interplane 
coupling of the multilattice compound Y2Ba4Cu70i5by 
modeling it with two inequivalent Hubbard layers cou- 
pled by an interlayer hopping t j_ . If the antifcrromagnetic 
correlation length is less than 1-2 lattice spacings, which 
happens for high temperatures or for a Fermi surface with 
suppressed nesting, we find that the two inequivalent lay- 
ers are disconnected and keep their individual magnetic 
properties. However, once the antiferromagnetic corre- 
lations in the layer with smaller charge carrier concen- 
tration is sufficiently large, the single-particle excitations 
for momenta close to the Fermi surface and, in partic- 
ular, around the momenta close to (tt, 0) as well as the 
magnetic excitations of the two layers are strongly con- 
nected. In this situation, the whole system reacts mag- 
netically as a single system, despite its inhomogeneous 
charge density. Due to the strong magnetic correlations 
in the underdoped layer the magnetic in-plane order in 
the nominally overdoped layer is stabilized via interplane 
magnetic coupling. Thus, the layer with lower charge 
carrier concentration acts like an external staggered field. 
For suitable Fermi surfaces and low enough temperatures, 
this mechanism makes the magnetic dynamics of both 
planes indistinguishable. By choosing parameters which 
lead to a moderate equalization effect, we have qualita- 
tively reproduced the salient features of the in-plane and 
out-of-plane spin-spin relaxation times observed in NMR 
measurements on Y2Ba4Cu70is, as compared with the 
ones of its constituent compounds. 
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